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PATENT 
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Multi-Frequency Pseudolites for Carrier-Based 
5 Differential-Position Determination 

The present invention relates generally to positioning systems 
for vehicles. It relates specifically to precisely determining positions for any 
land, sea, air, or space vehicle where existing satellite navigation methods 
such as GPS are inadequate due to limited visibility of its satellites. 

10 

This application claims the benefit of the following application: 
U.S. Patent Application No. 60/178,01 1, entitled, "GPS 

Performance in Deep Open Pit Mines using Pseudolites," filed January 24, 

2000, naming Clark E. Cohen et al. as inventors, with Attorney Docket No. 
1 5 P-68861 /AJT/LM and commonly assigned to IntegriNautics Corp. of Menlo 

Park, California. 

U.S. Patent Application No. 60/1 78,01 1 is incorporated by 

reference herein. 

20 Background 

Global positioning system (GPS) satellites transmit two coherent 
signals at different frequencies generated from a common time and 
frequency source. A receiver demodulates these signals to recover the 
underlying carrier frequency signal. The coherency of the two carriers 

25 facilitates the recovery of integer wavelength ambiguities in applications 
using carrier-based differential techniques. Many land applications of GPS 
(survey, ionospheric research, geophysical surface motion monitoring, etc.) 
use this technique. Typically, some residual ambiguity results from using the 
two coherent signals to resolve the integer wavelength ambiguities, and 

30 other methods are necessary to accelerate and/or confirm the accuracy 
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of the ambiguity resolution results. 

In a number of applications of carrier-based differential GPS 
(aircraft automatic landing and open-pit mine equipment positioning, for 
example), pseudolites augment the GPS satellite signals. In the landing 
5 application, the pseudolites facilitate the integer wavelength ambiguity 
resolution by incorporating the rapid geometry change due to the motion 
of the airplane close to the nearby pseudolites into the solution for the 
integer ambiguities. Cohen at aL U.S. Patent No. 5,572,218, describes such 
an application. 

10 In the open-pit mining application, for example, pseudolites 

augment the satellites obscured from view by high mine walls and cliffs. 
Single-frequency pseudolites used in this fashion provide additional 
code-phase measurements, but they do not aid in rapidly resolving carrier 
cycle ambiguities: The motion of the vehicles (shovels, trucks, crawlers, 

15 etc.) does not provide large geometry changes with respect to the 
pseudolites in a short period. 

Accordingly, an object of this invention is to rapidly resolve 
integer ambiguity — even without significant vehicle motion relative to the 
pseudolites. 

20 Another object of the invention is to enhance the integrity and 

speed of the pseudolite technique described above. 

Still another object of the invention is to leverage conventional 
GPS equipment, including GPS receivers and pseudolites, to reduce the 
cost of a system. 

25 These and other goals of the invention will be readily apparent 

to one of skill in the art on reading the background above and the 
description below. 
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Summary 

Herein are described apparatus and methods for rapidly 
resolving integer ambiguities in position determination. An embodiment of 

5 the invention includes a reference system, augmented with multi- 
frequency pseudolites using a carrier phase differential GPS 
implementation, and a mobile system. 

In one embodiment, the reference system may include one or 
more multi-frequency pseudolites, one or more multi-frequency reference 

1 0 receivers, a data link that is either stand-alone or built in to the pseudolites, 
and the associated antennae for each of these elements. The 
components of the reference system may be stationary. 

The mobile system may include a multi-frequency receiver and 
its associated antennae. Because the mobile systems may passively 

1 5 receive information, an unlimited number of mobile systems may be 
included in any given embodiment of the invention. 

A multi-frequency pseudolite uses a single frequency source to 
synthesize all of the multiple carrier frequencies and all of the multiple base 
band signals modulated onto carrier frequencies for transmission. The 

20 relative timing of all of the carriers and base band signals is constant and 
stable so no unknown timing drift exists between any of these signals. 

In a preferred embodiment, the modulation is a bi-phase 
pseudo random noise (PRN) sequence from the same family of Gold 
Codes used by the GPS satellites and the same code is used for all of the 

25 several carrier frequencies for any particular pseudolite. Different codes 
are used for the different pseudolites. The modulation timing is continuous 
for each pseudolite but the actual transmission of the modulated signals is 
pulsed to so that the pseudolites do not interfere with each other at the 
various system receivers. 

30 In the preferred embodiment, the multi-frequency pseudolite 

includes a multi-frequency reference receiver for synchronizing the 
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transmitted pseudoiite signals with the GPS system. This establishes 
coordination between all of the pseudolites in a given local 
implementation so that they can all be programmed to not overlap their 
pulsed transmissions. The integrated GPS receiver also surveys the location 

5 of the multi-frequency pseudolites. 

In the preferred embodiment wherein each pseudoiite is 
joined with a multi-frequency reference receiver, the receiver phase tracks 
all available GPS satellite signals and its associated pseudoiite signal and 
sends that data over the pseudolite's ranging signal. The mobile receiver 

10 then acquires this data directly from the RF ranging signals, similar to the 
way GPS sends satellite information at 50 baud as is known in the art. In 
other embodiments, where the multi-frequency reference receiver is 
independent of the pseudolites, the receiver not only phase tracks all 
available GPS satellites but all available multi-frequency pseudolites as 

1 5 well. The data in this case is relayed to the mobile receivers via a separate 
radio communication link. 

In one embodiment, the mobile receivers have functionally 
equivalent hardware and software as the reference receivers. They also 
receive through their associated antennae, to the extent possible, all of the 

20 same pseudoiite and GPS satellite signals received by the reference 
receivers but limited by line-of-sight restrictions due to surrounding 
obstacles. They also have software for assimilating the reference receiver 
phase information received over the data link from the reference receiver. 
The mobile receiver software combines the relative phase information for 

25 all of the pseudoiite and GPS satellite signals it has itself measured with the 
relative phase information for all of the pseudolites and GPS satellites that it 
has received over the data link from the reference receiver to then 
determine the wavelength integer ambiguities, thereby computing its 
position relative to the reference receiver position to an accuracy with 

30 errors on the order of only a fraction of a wave length at the LI (1575.42 
MHz) frequency. In most instances, the receiver may confirm, supplement 
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or substitute those calculations with additional calculations accurate to a 
fraction of a wavelength at each of the other multiple frequencies 
received at the mobile receiver. By using the multiple frequencies from the 
severa | pseudolites in view, the resolution of the wavelength integer 
5 ambiguities is nearly instantaneous (within one measurement cycle, 
typically 0.1 seconds or less) because it does not require a geometry 
change to accumulate due to either mobile receiver or GPS satellite 
motion to resolve the ambiguities. 

i o Brief Description of the Drawings 

Figure 1 illustrates an open-pit mine incorporating an 
embodiment of the invention. 

Figure 2 illustrates an open-pit mine incorporating another 
embodiment of the invention. 
1 5 Figure 3 illustrates an indoor positioning system incorporating 

an embodiment of the invention. 

Figure 4 illustrates a multi-frequency pseudolite according to 
one embodiment of the invention. 

Figure 5 illustrates a multi-frequency receiver according to one 
20 embodiment of the invention. 

Figure 6 shows a performance matrix for all combinations of 
satellites and pseudolites. 

Figure 7 is a table listing candidate frequency sets for three- 
and four-Frequency systems incorporating an embodiment of the invention. 
25 Figure 8 shows the cycle ambiguity probability for a single- 

frequency system as a nominal case. 

Figure 9 shows some benefits of the dual-frequency (or 
widelaning) GPS receiver. 

Figures 10 and 11 show that even with no motion in the system, 
30 the three and four-frequency designs can reduce the probability of 
selecting an incorrect cycle ambiguity. 
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Description of Specific Embodiments 

System Configuration 

Figure 1 illustrates an open-pit mine incorporating an 
5 embodiment of the invention. The mining application may include a pit 4, 
a mining vehicle 5, a global navigation satellite system 6 of the art and a 
multi-frequency pseudolite system 7. 

The pit 4 has walls 41 and a rim 42. 

The mining vehicle 5 includes a MFPS receiver 2. An MFPS 
10 receiver 2 recognizes and uses signals from the MFPS 7. Anything using a 
MFPS receiver 2 is a "user" of the MFPS 7, as the term is used herein. 
Typically, a user carries a multi-frequency mobile receiver 2 to determine its 
position. 

The multi-frequency pseudolite system 7 may include 

1 5 multi-frequency pseudolites 1 -1 :1 -N, where N is the number of pseudolites. 
The pseudolites may be to ensure that four or more pseudolites 1-i (i <= N) 
are available and provide good geometric dilution of precision (a factor of 
1 5 or less, preferably less than 5) to a user in the work space below — 
around the rim 42, for example. 

20 (Where no satellites are visible to a user, the invention may 

operate using only pseudolite signals, as described herein.) 

The satellites 3-l:3-K of any available global navigation 
satellite system (GNSS) may be incorporated. GPS is the most likely, but 
GLONASS may be included, as well as any future such systems. The 

25 satellites 3-1:3-K of the GNSS system 6 may respectively broadcast signals 
10-1:10-K and 11-1:11-K. The signals 10-1:10-K travel along a straight line 
towards a pseudolite 1-i of the pseudolite system 7 while the signals 11- 
1:1 1-K travel along a straight line towards a vehicle MFPS receiver 2. 

The surrounding terrain 41 may block some portion of the GNSS 

30 signals 11-1:11 -K, creating the need for the pseudolites. The mobile 
receiver 2 may use the remaining portion of the GPS signals 11 -1:1 1-K to 
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improve the ultimate position solution. 

Another portion of the GNSS signals 10-1:10-K may be used to 
survey the locations of the highest-elevation pseuclolites 1-j and to 
synchronize the system. Lower-elevation pseudolites (pseudolites 1-2 and 
5 1-3, for example) may need signals from upper-elevation pseudolites 
(pseudolites 1 -1 and 1 -N, for example) for survey and synchronization 
purposes. . 

A Multi-frequency Pseudolite 

1 0 Figure 4 illustrates a multi-frequency pseudolite 1 -i according 

to one embodiment of the invention. A multi-frequency pseudolite 1-i may 
include a receive antenna 217, a multi-frequency reference receiver 201, a 
control processor 202, a reference oscillator 203, a code generator 204, a 
pulse generator 205, multiple signal generators 210, a combiner 208 and a 

1 5 transmit antenna 207. These components may be interconnected as 
illustrated and as described herein. 

The signal generators 210 may produce GNSS-like radio- 
frequency (RF) signals for the pseudolite 1-i to broadcast. 

The signal generators 210 in the pseudolite 1-i may be any 

20 number. To operate in license-free RF spectrum bands currently available, 
the number may be four. 

Figure 7 is a table listing candidate frequency sets for three- 
and four-frequency systems incorporating an embodiment of the invention. 
For a four-frequency system, two frequencies may be in the 902.0 - 928.0 

25 MHz ISM band , and two frequencies may be in the 2400.0 MHz - 2483.5 
MHz ISM band. So as to not disrupt the normal use of any incorporated 
GNSS (GPS, for example), preferably none of the signal generators operates 
on any of the GNSS frequencies (LI (1575.42 MHz) and L2 (1227.60 MHz) for 
GPS, for example). 

30 A signal generator 210 may include a phase-lock loop (PLL) 

211, a filter 212, a microwave VCO 213, a mixer 214, a second filter 215 and 
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a switch 216. These components may be interconnected as illustrated and 
as described herein. 

The PLL 211 of a signal generator 210 may control the 
frequency produced by that signal generator 210 (the processor 202 
5 controlling the PLL 211). The reference oscillator 203 may drive the PLL 21 1 . 
The reference oscillator may be a low-cost temperature-controlled crystal 
oscillator (TCXO). The lower-frequency TCXO may stabilize the 
microwave-frequency carrier signal through the phase-lock loop feedback 
circuit 211, 212, 213. 

1 0 Where the GNSS is GPS, the code generator 204 may generate 

a gold-code (GPS C/A code) spreading pattern that the mixer 214 mixes 
with the carrier. Codes other than the GPS C/A codes may be used for the 
spreading function. Because of the economics that flow from using 
existing GPS chipsets, however, the GPS C/A codes or the GPS P-codes are 

1 5 the most likely choices for the spreading code. 

The code generator 204 may also modulate low-rate data (50 
to 1 k baud) on the carrier. The modulated data may include the surveyed 
location of the pseudolite, the pseudolite l-l's reference-oscillator offset 
from GPS time, and the satellite differential phase corrections and 

20 ephemeredes. Broadcasting all of this data on the pseudolite's ranging 
signal (the multiple RF signals) makes a separate differential reference 
station and subsequent communication link unnecessary. Multiple m 
separate data channels allows data to be sent in parallel (e.g., one byte 
can be distributed over four signal channels in two-bit chunks). This permits 

25 data rates up to m kbaud. Where m is four, the bandwidth is sufficient to 
accommodate all of the above data. 

The filter 215 filters the mixed carrier, code, and data, and 
under the control of the pulse generator 205, the switching circuit 216 
passes the filtered mix. Pulsing a pseudolite signal at a short duty cycle 

30 (say, 5-10%) and a relatively high power (say, 1-10 milliwatts) may provide 
benefits when a pseudolite transmits on a GNSS frequency. For example, 
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pulsing may enable a pseudolite to operate without significant impact on 
GNSS satellite signals. Also, pulsing may increase the useful spatial range of 
the pseudolite signal. Without pulsing, the range-ratio over which the 
pseudolite can be used is roughly 3:1 ; with pulsing, the range-ratio is 
5 greater than 1000:1. 

(Where the multiple signals broadcast by the multi-frequency 
pseudolite are not on the GNSS frequencies, pulsing may not provide the 
first benefit — but in this instance, the first benefit is no longer required.) 

The combiner 208 may sum the multiple signals from the signal 

10 generators 210, and the transmit antenna 207 may then broadcast them. 
The antenna 207 may be a wide-band antenna or may be separate 
narrow-band antennae. Whatever the antenna arrangement, coincident 
phase centers are preferred for all the frequencies involved. Otherwise, the 
position solution may be more challenging to calculate than presented 

15 herein. 

A pseudolite #1-1 may include a multi-frequency reference 
receiver 2 capable of positioning from multi-frequency pseudolites. First, 
such a receiver 2 enables the pseudolite 1 -i to automatically survey its 
location and then broadcast that information via the data message. 

20 Second, the receiver 2 enables the pseudolite 1 -i to collect satellite 

differential correction data and ephemeredes (sent via the data message). 
Third, the receiver 2 provides synchronization for triggering the pulse 
generator 205. Since the receiver 2 can use other multi-frequency 
pseudolites 1-i to determine its position, the pseudolite 1-i may be placed 

25 in an area of insufficient satellite coverage yet still be fully operational. 

For indoor applications of the invention, the integrated receiver 
201 may be omitted by connecting all the pseudolites to a common 
oscillator. The locations of the pseudolites in an indoor setting are surveyed 
by means other than GNSS, and the survey data entered — manually, for 

30 example — into the pseudolites for re-broadcast to users. 
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A Multi-frequency Receiver 

Figure 5 illustrates a multi-frequency receiver 2 or 201 
according to one embodiment of the invention. The receiver 2 of Figure 5 
5 is the counterpart to the pseudolite 1-i of Figure 4 , and both assume that 
GPS is the operative GNSS. 

(It is worth noting here that a multi-frequency receiver in the 
context of this patent pertains to a receiver capable of interpreting the 
signals of a multi-frequency pseudolite as described herein. A multi- 
10 frequency receiver is different from a "dual-frequency" or "L1/L2" GPS 
receiver, which is a well-known technology.) 

The receiver 2 may include a receive antenna 317, a splitter 
316, multiple up/down converters 315, respective multiple GNSS receivers 
301 and a navigation processor 305. These components may be 
15 interconnected as illustrated and as described herein. 

A standard LI GPS receiver with a frequency converter and 
modified internal software can interpret a non-U frequency signal just as if 
it were broadcast on LI . Thus, a multi-frequency receiver may be built up 
from a set of such LI receivers and frequency converters. Such a receiver 
20 may include five conventional GPS receivers 305 under control of the 

navigation processor 305. Each of four of the receivers 305 may handle a 
respective one of the four signals transmitted by the pseudolites 1-i . A 
frequency converter 315 between the antenna 317 and the front end of a 
GPS receiver 301 may modulate the incoming signal up or down to GPS LI 
25 so the receiver 301 can work with the signal as if it were broadcast on LI . 
The fifth GPS receiver 301 may directly measure the satellite signals and 
does not require a frequency converter. 

More generally, the multi-frequency receiver 2 may include 
m + 1 conventional GNSS receivers 305 under control of the navigation 
30 processor 305. Each one of m of the GNSS receivers 305 handles a 
respective one of the m signals transmitted by the pseudolites 1-i . The 
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frequency converter 315 modulates the incoming signal up or down to a 
predetermined standard GNSS frequency so the standard receiver 301 can 
work with the signal as if it were broadcast on that predetermined standard 
GNSS frequency. 

5 Like the pseudoiite transmit antenna 207, the receiver antenna 

317 may be a wide-band antenna or separate narrow-band antennae, in 
any event with a coincident phase center for all frequencies involved. A 
low-noise amplifier (LNA, not shown) may amplify the received signal and 
the five-way (m + 1 ) splitter may supply a signal to each of the GPS 

1 0 receivers 301 . 

A GPS receiver 301 may include a RF front-end and sampler 
302, correlator channels 303 and a phase processor 304. 

The GPS receivers 301 may not process a position solution as a 
conventional GPS receiver does but may collect phase measurements in 

1 5 the form of code and carrier phases. The RF front-end and sampler 302, 
correlator channels 303, and phase processor 304 perform this function as is 
well known in the art. The navigation processor 305 may run a process 307 
that collects code and carrier phase measurements 306 from the receivers 
and sends them to the navigation solution process 308. The receiver 2 

20 outputs the navigation solution. 

In processes 311 and 314, the navigation processor 305 may 
use the navigation solution to command the GPS receivers as to which PRN 
codes to search for and which Doppler offsets to use in searching for 
additional satellites and pseudolites. 

25 

Additional Embodiments 

Figure 2 illustrates an open-pit mine incorporating another 
embodiment of the invention. The embodiment may include a pit 4, a 
mining vehicle 5, a global navigation satellite system 6 of the art and a 
30 multi-frequency pseudoiite system 7, as well as a multi-frequency reference 
station 14 separate from the multi-frequency pseudolites of the MFPS 7. 
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The reference station 14 includes a reference receiver 2. (This embodiment 
may prove technically easier to implement in a real application of the 
invention than that of Figure 1). 

With line-of-sight visibility to all of the same pseudolites and 
5 satellites that all of the mobile receivers have available, the reference 
station 14 may serve the entire work space. The reference receiver 2 may 
collect differential code and carrier phase information and satellite 
ephemeredes and distribute these to users' mobile receivers 2 via a radio 
communications link 13 independent of the pseudolites. 
1 0 Figure 3 illustrates an indoor positioning system incorporating 

an embodiment of the invention. In this instance, no satellites are visible. 
An oscillator 15 connected through one or more cables 16 to pseudolites 
1-i may synchronize the pseudolites 1-i. 

1 5 Position Solution Using Multi-frequency Pseudolites and Receivers 

A multi-frequency GNSS receiver 2 can calculate its position by 
measuring the code and carrier phases transmitted by a set 7 of 
multi-frequency pseudolites 1-i and processing these phase measurements 
with the algorithms described herein. The position calculation differs from 
20 the conventional GNSS position determination. (While any GNSS may be 
used to illustrate position solving, the following explanation assumes the 
GPS.) 

Existing GPS receiver technology typically provides 
code-phase measurements with meter-level noise figures and carrier-phase 
25 measurements with centimeter-level noise figures. With four or more 
multi-frequency pseudolites in view and with good geometry, the 
multi-frequency GPS receiver can employ carrier-phase measurements 
and immediately achieve centimeter-level position solutions. 

In embodiments described above, each pseudolite 1-i 
30 includes a reference receiver 2. An included reference receiver 2 can be 
assumed to be zero range from its respective pseudolite 1-i and to provide 
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a direct measure of the difference between the local oscillator 203 and 
GPS time. From these assumptions, it follows that all phases between the 
pseudolite and a user can be considered differentially corrected and the 
only time bias in the system is between the mobile receiver and GPS time. 
5 The parameters of interest in the position calculations are as 

follow: 

R ( . True range between the mobile user and pseudolite /. 
p { Precise range derived from phase measurements. 
10 f- } Frequency of carrier signal /,/ E {1 ...N f }, 

A } Wavelength of carrier signal j, A } = c/f^, where c is the speed of 
light. 

</\ Code phase to pseudolite / for frequency j (differentially 
corrected). 

15 ifr X] Carrier phase to pseudolite / for frequency / (differentially 

corrected). 

T Range-equivalent of the offset between the mobile-receiver 
clock and GPS time. 

M Vj Code-phase cycle ambiguity from pseudolite / to mobile receiver 
20 for frequency j. 

N l} Carrier-phase cycle ambiguity from pseudolite / to mobile 
receiver for frequency /. 

N f Number of frequencies implemented in the ranging signal. 

C7 cd Variance of the code phase measurements (assumed same 
25 for all carrier frequencies) . 

a } Variance of the carrier phase measurements for frequency /. 

w cd Code phase noise as zero-mean Gaussian with variance a cd . 

Wj Carrier phase noise as zero-mean Gaussian with variance <3 r 

30 The first step of the positioning algorithm determines the precise 

ranges between the mobile user and each of the pseudolites. "Precise 
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range" in this instance means that the range accuracy is established to 
better than one wavelength of the highest frequency employed. The 
second step extracts the position coordinates from the precise-range 
measurements. The second step is fairly well understood in the art — 
5 employing a non-linear least-squares algorithm, for example — and so only 
the firsf step is covered in detail. 

The range R } from the pseudolite to receiver station is the 
magnitude of the difference between vectors P and Q, and is related to 
the code phase measurement ^ by equation (1 ): 

10 

As noted above, for the configuration under consideration, is a 
differential code phase, and 7ls the offset between the mobile receiver's 

1 5 clock and GPS "true time." M }j is the code cycle ambiguity from the 
pseudolite / to the mobile receiver for frequency j. 

For most applications of the present invention, M, can safely be 
assumed to be zero because the scale of the invention will be less than 
one code phase cycle (300km) in any dimension. This simplifies the code 

20 measurement equation to equation (2): 

(2) 

The carrier-phase measurements are similar in form, but the cycle 
ambiguities N ti must be resolved. Again, due to the co-location of the 
25 reference receiver with the pseudolite, these are differential carrier phases 
and the time bias is between the mobile receiver and true time: 
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yr v =R+r + A,N 9 + Wj 



(3) 



The time bias, T, is common to all measurements and can be omitted from 
the equations for the time being. It will be re-introduced after the precise 
range is resolved: 



'a (4) 



10 

Resolving the precise ranges involves superimposing all the 
probability density functions for the various carrier cycle ambiguities from a 
given pseudolite, centered around the code phase measurement and 
enveloped by the code-phase probability density function. The probability 

1 5 density function for a given carrier can be viewed as a comb with spacing 
between the "teeth" equal to the wavelength of that carrier. By 
superimposing combs of different spacing (the probability density functions 
for the various carrier signals), only one set of "teeth" overlap, representing 
the most probable solution. The width of the comb (and hence the total 

20 number of "teeth") is bounded by the probability density function for the 
code phase measurement. 

The starting point for the resolution process is the approximate 
range provided by the code phases. If there are N r frequencies involved, 
then there will be N f code phase measurements from a given pseudolite. 

25 <p l0 is the average of the code phases between pseudolite / and the mobile 
receiver: 
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(6) 



Since the raw carrier phase measurements contain an arbitrary 
cycle ambiguity, this ambiguity is stripped off from the fractional phase and 
replaced with the closest number of integer wavelengths as provided by 
5 the code phase estimate $ 0 . The derived quantity is the centered carrier 
phase / dj . In the following equation, the modulus operator strips off the 
integers from the raw carrier phase measurement to leave a fractional 
phase, and the floor function appends the range of the closest integer to 
the code phase estimate <p ]Q \ 

10 

The centered carrier phase / cjj serves as the location to center 
the probability density function for carrier /. After all of the probability 
1 5 density functions are superimposed, the highest probability integer 

indicates an offset from the original code phase estimate (f> [Q . This offset is 
applied to the original code phase estimate to achieve the precise range 

Pr 

Sampled sequences of the probability density functions for all 
20 the carriers can be generated and then multiplied together to determine 
the highest probable cycle ambiguity solution. The range (length) of the 
sampled sequence is determined by the code phase variance <J cd . The 
sequence is preferably at least 2-3 variances wide. The resolution of the 
samples is preferably finer than the resolution of the highest frequency 
25 carrier phase measurement (typically less than a centimeter). The sampled 
probability density function for pseudolite /, frequency / is represented as 
p iJk , where k is the sample index and is directly related to the range by 
k *dr, where dr is the sample resolution. Assuming the carrier phase 
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measurement noise is zero-mean Gaussian with variance q, the probability 
density function p Jk is the superposition of independent Gaussian 
distributions, each with variance and separation between centers of one 
carrier wavelength A r Additionally, the entire distribution is shifted by the 
5 difference between the centered carrier phase, / cij , and the code phase 
estimate $ 0 . p Jk can be represented mathematically as equation (8): 



{k-dr-A-j-n+Pcj-fa,) 2 



2o] 



(8) 



1 0 Superimposing the probability density functions for all of the carriers j for a 

given pseudolite ;', and enveloping the code probability density function 

results in the probability density product sequence, P lk . 

P jk is characterized by a unique maximum value or spike that 

indicates the location in the sequence of the correct integer solution. P ik is 
1 5 formed by multiplying the code probability density sequence and the p ljk 

sequences together: 



(bdrf N f 

P*=e Tp ijk ^ 



20 By finding the index, m, of the maximum value of P !k , the offset between the 
code phase estimate and the precise range is determined to be m *dr. 
The precise range is thus computed as in equation 10: 



Pr0- m ' dr 



(10) 



25 



The precise ranges between all pseudolites and the mobile 
receiver can be found in this manner. It should be noted that the precise 
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range established here is actually the precise range plus the time bias that 
was omitted early in the derivation. The procedure for extracting the 
position solution and the time bias from a set of precise range 
measurements such as these is well known in the art. 
5 For universal application of the invention, preferably at least 

five pseudolites are in view at all times. This arrangement provides the 
highest performance available from the system. It is useful to know, 
however, the expected performance when fewer than five pseudolites are 
in view. 

} o Figure 6 shows a performance matrix for all combinations of 

satellites and pseudolites. When at least five pseudolites are in view, the 
user receiver can expect instantaneous integer acquisition, full integrity, 
and the highest level of accuracy. If four pseudolites and at least one 
satellite are in view, the system has similar performance — except that the 

1 5 integrity is not guaranteed without resorting to motion-based algorithms to 
independently resolve the cycle ambiguities for the satellites. These 
algorithms may take 10 to 20 minutes. 

The need for integrity is highest when the navigation sensor is 
used in feedback control of safety-critical machinery. For existing man-in- 

20 the-loop operations this is not a firm requirement, but it can serve to 

improve the safety of these operations nonetheless. Full integrity monitoring 
does provide the means for achieving higher levels of automation in many 
systems. Situations in which there are fewer than four pseudolites but at 
least five ranging sources in total can achieve full accuracy and integrity, 

25 but resort to satellite motion to resolve cycle ambiguities. 

The expected performance of any set of frequencies to 
uniquely discern the correct integers can be determined by computing the 
probability of a cycle ambiguity occurring over the range of the code 
phase accuracy. An example of a candidate frequency set for each of a 

30 three-frequency and a four-frequency architecture are listed in Figure 7. 
The table lists a fourth band designated for ISM in the 5725 - 5875 MHz 
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range which remains as an option if the other bands are overcrowded by 
other radio equipment in the area of the application. 

The ability to reject erroneous integer cycle solutions using 
these frequency sets can be compared with the conventional GPS 
5 widelane technique by observing the graphs of the cycle ambiguity 
probability versus range. The cycle ambiguity probability for a single- 
frequency system is shown in Figure 8 as a nominal case. 

This semi-logarithmic graph shows that carrier phase alone 
does not distinguish between being at range 0 as opposed to range 0.19 

10 m, 0.38 m, 0.57 m, etc. Each spike (except for the one at zero) is a 

potentially wrong selection for the carrier cycle resolution algorithm. The 
only discernment of range comes from the code phase probability 
envelope, which in this case is modeled conservatively as a normal 
distribution with a standard deviation of 3m. The width of each spike 

15 represents the carrier phase noise, modeled here as a normal distribution 
with a standard deviation of 4% of a carrier cycle, or 0.75 cm for LI . The - 
ability to suppress cycle ambiguities in a multi-frequency system is a 
function of code and carrier phase noises, as well as the frequencies 
themselves. 

20 Figure 9 shows some benefits of the dual-frequency (or 

widelaning) GPS receiver. The superimposed cycle ambiguity probability 
functions of LI and L2 reduce the possible integer solutions substantially, but 
not enough to completely obviate the need for satellite motion to discern 
the correct cycle. This is therefore still not adequate for a high-integrity 

25 motionless solution. 

Even with no motion in the system, the three and four- 
frequency designs can reduce the probability of selecting an incorrect 
cycle ambiguity to less than 1 in 10 and 1 in 500, respectively, as shown in 
Figures 10 and 11. This substantial improvement over the dual-frequency 

30 system provides the ability to resolve position to better than 10 cm in one 
measurement period. The 1 in 10 possibility for a cycle ambiguity error in 
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the three-frequency system is deemed marginal performance, although it is 
worth noting that redundancy provided by five or more signal sources (i.e. 
any combination of satellites and pseudolites) establishes the ultimate 
integrity of the navigation solution. Integrity monitoring via redundant 
5 measurements readily discerns the 1 in 500 event of an incorrect cycle 
ambiguity in the four-frequency case. 

The present invention combines and extends the pseudolite 
technique with the dual frequency satellite technique to create apparatus 

1 0 and methods for rapid integer ambiguity resolution. The invention further 
enhances the integrity and speed of the technique by increasing the 
number of coherent signal sources to three, four or even more frequencies. 
The invention uses multi-frequency pseudolites to accelerate and facilitate 
the resolution of integer wavelength ambiguities in carrier-based differential 

15 GPS applications. Consequently, integer ambiguities can be determined 
extremely quickly and without any residual uncertainty as to the accuracy 
or fidelity of the wavelength ambiguity determination. This eliminates 
augmenting the resolution process with other techniques requiring satellite 
or vehicle motion. (Satellite motion takes time to produce a useful 

20 geometry change, while vehicle motion may not be useful or practicable 
in applications using slow or ponderous equipment.) 

In one embodiment, the invention is a multi-frequency 
pseudolite that provides the signal-in-space needed by a receiver to 
determine the cycle ambiguity of the microwave carrier frequencies. The 

25 cycle ambiguity is determined with very high and quantifiable accuracy 
and integrity. 

In another embodiment, the invention is a 
microwave-frequency receiver that leverages a multi-frequency pseudolite 
signal, as well as standard GPS satellite signals, to rapidly determine its 

30 position to the accuracy typical of carrier-phase differential GPS 
techniques. 
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In another embodiment the invention leverages conventional 
GPS equipment, including GPS receivers and pseudolites, to reduce the 
cost of the overall system. 

5 Indeed, the invention now being fully described, many 

changes and modifications that can be made thereto without departing 
from the spirit or scope of the appended claims will be apparent to one of 
ordinary skill in the art. 

10 
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